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ABSTRACT 


The steppe viper was rediscovered in Georgia after 75 years. A comprehensive analysis of external morphology, 
altitude gradient in habitats, typology of biotopes and genetic analysis revealed a high degree of similarity of popula- 
tions of the steppe vipers from Azerbaijan, known as Pelias shemakhensis Tuniyev et al., 2013, and from East Georgia. 
These data were used for comparative study and description of a new subspecies — P. shemakhensis kakhetiensis ssp. 
nov. The subspecies name is after historical region of Georgia — Kakheti, where a large part of the range is located. 
In the pattern of recent distribution of P. shemakhensis, there are common habitat and climatic characteristics in 
Georgian and Azerbaijan parts of this its range. Its position among the species complex and relations with other taxa 
of this complex are discussed. Based on the results of the cluster and discriminant analyses, P. eriwanensis and P. lo- 
tievi should be given a subspecies rank, whereas P. shemakhensis clearly deserves a species rank. Results of the genetic 
analysis are opposite: P. shemakhensis and P. eriwanensis are considered as sister species. We continue to share an idea 
about autochthonous origin for small shield-head vipers of the Caucasus, taking into account their ancient origin in 
the Caucasian Ecoregion and astonished variety of forms of the “kaznakovi”-complex and “ursinii”-complex both on 
the Great Caucasus and in the Transcaucasia. 
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PE3IOME 


Спустя 75 лет степная гадюка вновь была обнаружена Ha территории Грузии. Комплексный анализ внеш- 
ней морфологии, высотного градиента обитания, типологии биотопов и генетический анализ показали 
высокую степень сходства популяций степных гадюк из Азербайджана, известных как Pelias shemakhensis 
Tuniyev et al., 2013, и Восточной Грузии. Полученные материалы были использованы для сравнительного 
анализа и описания нового подвида шемахинской гадюки — Р. shemakhensis kakhetiensis ssp. nov. Подвид на- 
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зван по географическому наименованию исторической области Грузии — Кахетии, в которой расположена 
большая часть ареала. Рассматривая современное распространение P. shemakhensis, обращает внимание сход- 
ство грузинского и азербайджанского кластеров ареала по биотопическому и климатическому факторам. 
Рассмотрено положение нового подвида в комплексе видов степных гадюк и обсуждаются его отличия от 
других таксонов этого комплекса. По результатам кластерного и дискриминантного анализов, Р. eriwanensis 
и P. юйе выступили скорее B подвидовом ранге, тогда, как P. shemakhensis четко отделяется в видовом ранге. 
По результатам генетического анализа получен противоположный результат — Р. shemakhensis и Р. eriwanensis 
являются сестринскими видами. Мы продолжаем придерживаться мнения 00 автохтонности мелких 
щиткоголовых гадюк Кавказа, учитывая их древность в Кавказском экорегионе и поразительное разнообра- 


aue форм "kaznakovi" -complex и “ursinii”-complex как на Большом Кавказе, так и в Закавказье. 


Ключевые слова: Грузия, новый подвид, Pelias shemakhensis, степная гадюка 





INTRODUCTION 


A guestion about the occurrence of the steppe 
viper in Georgia remained open from the middle of 
the last century. Previously, this adder seems to be 
locally abundant there: Dzhanashvili (1951) noted 
that in 1944 13 specimens were killed in the distance 
of 15 km of the road in the Shirak Steppe. Later, 
nobody has confirmed the occurrence of the species 
within Georgia. All special surveys of the vipers in 
Eastern Georgia had no successful results. For this 
reason, it was suggested that the steppe viper really 
inhabited Eastern Georgia earlier, but subseguently 
became extinct (Muskhelishvili 1968). Thus, the 
clarification of the possible recent occurrences of this 
species is actual (Tuniyev et al. 2013). 

In this regard, findings of the steppe vipers in 
2015 in a number of closely-spaced localities on the 
Shirak Plateau (Kakhetia, East Georgia) made by 
G.N. Iremashvili were of great interest, also in the 
aspect of their taxonomic position and relations with 
the recently described Pelias shemakhensis Tuniyev et 
al., 2013 (Tuniyev et al. 2013). A hypothesis about 
possible belonging of the vipers from Eastern Georgia 
to the same taxon was considered. 

Institutional and other abbreviations. Az — 
Azerbaijan; Ge — Georgia; MNKNU - Museum of 
Nature at V. N. Karazin Kharkiv National University, 
Kharkiv, Ukraine; SMG — State Museum of Georgia, 
Tbilisi, Georgia; SNP — Herpetological collection of 
the Sochi National Park, Sochi Russia; ZIRA — Sci- 
entific Center of Zoology and Hydroecology of the 
National Academy of Sciences of Republic Armenia, 
Yerevan, Armenia; ZISP — Zoological Institute of 
the Russian Academy of Sciences, Saint Petersburg, 
Russia. 


MATERIAL AND METHODS 


A total of 65 specimens of vipers from the Cauca- 
sian Ecoregion, related to the “ursinii”-complex were 
examined including seven specimens of the steppe 
viper from East Georgia. Data on morphology of the 
type specimen of Pelias shemakhensis (ZISP 21720), 
as well as ten specimens from Shemakha and two 
specimens from Demirchi (Azerbaijan), published by 
Kukushkin et al. (2012) were also used. 

77 specimens of vipers from East Georgia, Azer- 
baijan, Russia (Dagestan) and Armenia were used 
in statistical, cluster and discriminant analyses. 
Material is stored in the herpetological collections 
of the SNP, ZIRA and ZISP (see Table 1). Material 
is combined according to geographical localities into 
four samples: 1) Pelias eriwanensis (Reuss, 1933), Ar- 
menia; 2) Pelias shemakhensis, Shemakha, Demirchi, 
Azerbaijan; 3) Pelias cf. shemakhensis, Shirak Plateau, 
East Georgia; 4) East Pelias lotievi (Nilson et al., 
1995), Dagestan, Russia. 

Morphological analyses. The methods of tradi- 
tional morphology were used for morphological char- 
acters offered by Nilson and Andrén (2001) with some 
our modification (Table 2). To eliminate influence of 
age variation, comparison of adult and juvenile males 
and females were made separately only on meristic 
characters (Table 2). Material was studied statisti- 
cally using standard methods of variation statistics 
(Lakin 1990) and one of the methods of multidimen- 
sional statistics — Canonical Discriminate Analysis 
(CDA), and cluster analysis (Tyurin et al. 2003) by 
the package of STATISTICA 6.0 for Windows. 

DNA extraction, amplification and sequencing. 
23 specimens sampled from 14 localities (Fig. 1) were 
analyzed. Five sequences of Pelias ebneri (Knoep- 
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Table 1. Examined specimens of Pelias spp. See text for abbreviations. 
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Tree Genbank | . 
No. No. accession No. Species n Sample Locality Date 


ZISP Vicinity of Shemakha, 


21720 - - P. shemakhensis 1 2 Azerbaijani 14.04.1974 
SNP P. shemakhensis Khornabudji Fortress, Shirak 
1059 420 MG729486 kakhetiensis ssp. nov. 1 3 Plateau, Kakhetia, East Georgia 06.2016 
А Vicinity of the Chinkani Village, 
SNP 399 МС729487 Р shemakhensis ; | 
1052 402 MG729488 kakhetiensis ssp. nov. 3 pee Plateau, Kakhetia, East 03:2019 
eorgia 
; Vicinity of the Chinkani Village, 
SNP од.  mMmG729489 7 Shemakhensis 3 Shirak Plateau, Kakhetia, East 04.2016 
1053 kakhetiensis ssp. nov. Georgi 
eorgia 
ZIRA Tapchan Elakh, above Kakava- 
40715 - - P. eriwanensis 1 1  berd Castle, *Khosrov Forest" 06.2011 
Reserve, Ararat Marz, Armenia 
no P - = Р, eriwanensis 1 1 The same locality 06.2011 
ZIRA : : . 
40713 - - P. eriwanensis 1 1  Thesamelocality 06.2011 
ZIRA : : . 
40718 - - P. eriwanensis 1 1  Thesamelocality 06.2011 
p "n - - P. eriwanensis 1 1  Thesamelocality 06.2011 
ZIRA Jabrail, vicinity of Getatakh and 
- - P. eriwanensis 1 1 Lor villages, Sisian, Sjunik Marz, 12.05.2006 
40716 i 
Armenia 
ZIRA : : . 
40721 - - P. eriwanensis 1 1  Thesamelocality 12.05.2006 
ZIRA : : к 
40719 - - P. eriwanensis 1 1  Thesamelocality 12.05.2006 
ZIRA : : : 
40720 - - P. eriwanensis 1 1  Thesamelocality 12.05.2006 
s - - P. eriwanensis 1 1  Thesamelocality 12.05.2006 
SNP Я , Vicinity of Drakhtik Village, 
890 43 MG543319 P. eriwanensis 2 1 Areguni Ridge, Armenia 
Branch of Katar Mountain, 
SNP ; | above Verin Giratakh Village, 
894 42 MG543320 P. eriwanensis 1 1 Khustup-Katarsky Ridge, 11.05.2012 
Sjunik Marz, Armenia 
Uu 219 | MG543321 Р егіоапепѕіѕ 3 1  Thesamelocality - 
Near Asni Village, southern slope 
ZIRA | А Е d 8 and 
708 - - P. eriwanensis 2 1 of Urtz Ridge, Ararat Marz, 18.05.1939 
Armenia 
ABE - Reriwanensis gw V Ridge Ararat Marz, 10.05.1950 


812 Armenia 


Collector 


L.A.Erukh 


N. Seturidze 


GN. Iremashvili 


GN. Iremashvili 


A. Malkhasjan 


A. Malkhasjan 
A. Malkhasjan 
A. Malkhasjan 


A. Malkhasjan 


A. Malkhasjan 


A. Malkhasjan 
A. Malkhasjan 
A. Malkhasjan 


A. Malkhasjan 


19—20.05.2012 S. Tuniyev 


B. Tuniyev 


A. Malkhasjan 


S.K. Dal 


P.P. Gambarjan 
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Table 1. Continued. 
No Ts Genbank Species n Sample Locality Date Collector 
` No. accession No. 
ZIRA Sukhoy Fontan Settlement, 
615 - - P. eriwanensis 1 1  Akhtynsky (Razdan) District, 03.06.1936 — Izmailov 
Armenia 
SNP 170 | MG543322 
797 171 MG543323 P. eriwanensis 3 1 Arai-Ler Mountain, Armenia 08.2008 S. Balajan 
172 MG543324 
SNP А А bakh AL Р 
1049 45 MG543325 Р. eriwanensis 1 1  Nagorno-Karaba' 07.2012 .L. Agasjan 
ZIRA : : бе ae : : 
n/n 220 MG543326 P. eriwanensis 1 1 Sisian, Sjunik Marz, Armenia ? A.L. Agasjan 
Vicinity of Shari Village, : 
SNP 160 MG729490 P lotievi 3 4 Djufudag Ridge, Agulsky 21.08.2007 B S Tuniyev, 
733 Smee 4 S. B. Tuniyev 
District, Dagestan, Russia 
SNP aes : B. S. Tuniyev, 
738 - - P. lotievi 4 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP i a : B. S. Tuniyev, 
739 - - P. lotievi 5 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP - : B. S. Tuniyev, 
740 - - P. lotievi 1 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP "m : B. S. Tuniyev, 
743 - - P. lotievi 1 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP i os : В. S. Tuniyev, 
744 - - P. lotievi 1 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP а : B. S. Tuniyev, 
745 - - P. lotievi 1 4 The same locality 21.08. 2007 S. В. Tuniyev 
SNP = 5 B. S. Tuniyev, 
746 - - P. lotievi 1 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP bai ; B. S. Tuniyev, 
750 - - P. lotievi 1 4 The same locality 21.08. 2007 S. B. Tuniyev 
SNP Vicinity of Burshag Village, 
- - P. lotievi 1 4  Djufudag Ridge, Agulsky 21.08.2007  B.S. Tuniyev 
749 dn : 
istrict, Dagestan, Russia 
Vicinity of Ritlyab Village, : 
SNP - - P. lotievi 20 4  Charodinskiy Reservation, 19.08. 2007 B.S. Tuniyev, 
735 . S. B. Tuniyev 
Dagestan, Russia 
SNP ах Vicinity of Tidip Village, Shamil р 
764 - - P. lotievi 3 4 District, Dagestan, Russia 22.06.2008 Z. Ismailova 
SNP o Verkhnee Gakvari Village, А 
762 144 MG729491 Р. lotievi 3 =. Dagestan, Russia 02.07.2008 В. S. Tuniyev 
SNP un Pass Gumbashi, . 
1075 422 | MG729492 P lotieoi 1 = Karachay-Cherkessia Russia 23.08.2016 Kh. U. Aliev 
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Table 1. Continued. 








Tree Genbank | , 

No. Nu. bean No. Species Locality Date Collector 
SNP 243  MG729493 P lotievi Tareim Hollow, Erzi Reserve, 18.06.2014 К. Yu. Lotiev 
954 Ingushetia, Russia 

Vicinity of Saratovskaya : : 
SNP 492 MG729494 P renardi Settlement, Krasnodar Territory, 20.04.2012 WN: Timukhin, 
901 А S. В. Tuniyev 
Russia 
185 МС729495 Vicinity of Saratovskaya . 
S 186 MG729496 P.renardi Settlement, Krasnodar Territory, 06102012 85 Tuniyev, 
187 MG729497 Russia UU y 
Kirista Peninsular, northern 
SNP 201 MG543327 ; bank of Manych Lake, protected S.B. Tuniyev, 
782 204 MG543328 Prenardi area of Chernye Zemli Reserve, 1805-2009 S Taniyev 
Kalmykia, Russia 
SNP Я Kaban Mountain, vicinity of K.Yu. Lotiev, 
785 208 МС729498 P. renardi Kislovodsk Town, Russia 23.05.2009 S.B. Tuniyev 








Fig. 1. Localities of Pelias spp. tissue samples for DNA analyses: 1 — Arailer Mountain, Armenia; 2 — Sevan Lake, Drakhtik Village, 
Armenia; 3 — Sisian, Armenia; 4 — Vicinity of the Verin Giratakh Village, Bargushat Ridge, Armenia; 5 — Nagorno-Karabakh; 6 — Khorn- 
abudji Fortress, Georgia; 7 — Chinkani Village, Georgia; 8 — Shari Village, Dagestan, Russia; 9 — Verkhnee Gakvari Village, Dagestan, 
Russia; 10 — Erzi State Reserve, Ingushetia, Russia; 11 — Kaban Mountain, Stavropol Territory, Russia; 12 — Gumbashi Ridge, Karachay- 
Cherkessia, Russia; 13 – Manych-Gudilo Lake, Kalmykia, Russia; 14 — Saratovskaya Settlement, Krasnodar Territory, Russia. 
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Table 2. Abbreviations for morphological characters and indexes of vipers. 








No. Abbreviation Description 
1 Т1. Total length 
2 SVL Snout-vent length 
3 L.cd. Distance from cloaca to tail tip 
4 Pr. Number of preventral shields 
5 Ven. Number of ventral shields 
6 S.c. Number of subcaudal shields 
7 Ap. Number of apical shields 
8 R. Rostral index: ratio of height to width 
9 Pil. Distance from tip of snout to posterior margin of parietal shields 
10 Crown scales (C.s.) Number of shields, limited by frontals, canthals and supraoculars 
11 In Upper preoculars in contact with nasal, left/right (+/-) 
12 Can. Number of canthal shields 
13 Sa.1 Number of dorsal scales around neck 
14 Sg.2 Number of dorsal scales around midbody 
15 Sa.3 Number of dorsal scales around posterior part of body 
16 Supralab. Number of supralabial shields 
17 Sublab. Number of sublabials shields 
18 Ec. Number of shields around eye (left/right) 
19 ZZ Number of zigzag windings (left/right) 
20 Lor. Number of loreal shields (left/right) 
21 Lc. Head length 
22 Lt.c Biggest width of head 
23 Al.c. Height of head in parietal area 
24 Par. Parietals (hel/delad) 
25 Fr. Frontal (hel/delad) 
26 Nas. Nasal (hel/delad) 





pfler et Sochurek, 1955) were taken from GenBank 
(KC176745, FR727093-FR727096). Genomic DNA 
was isolated from tissues fixed with 9696 ethanol 
using standard salt extraction protocol using a lysis 
buffer and proteinase K, deproteinized with NaCl and 
precipitated with 96?6 ethanol (Miller et al. 1988). A 
segment (1154 bp) of the cytochrome b (cytb) gene 
was amplified using primers Cytb F1 and Cytb RC, 
additional internal primer Cytb F8 was used for se- 
quencing (Stumpel 2012). 


PCR conditions were as published in Zinenko et 
al. (2015). Each PCR included a negative control. 
The PCR products were purified on columns of an 
Omnix kit (Omnix, Saint Petersburg, Russia) and 
were sequenced in both directions using the BigDye 
Terminator Cycle Sequencing Ready Reaction Kit 
on an ABI PRISM 3130 (Applied Biosystems Inc., 
Foster City, CA, USA). Sequences were edited and 
assembled using program Sequencher 4.6 (Gene 
Codes Corporation, Ann Arbor, MI, USA http:// 
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www.genecodes.com) and aligned with CLUSTALW 
algorithm (Thompson et al. 1994) implemented in 
BioEdit (Hall 1999). 

Phylogenetic analyses. Tree reconstructions 
were performed using 28 specimens of Pelias spp. 
and 3 specimens as outgroup (Bitis arietans (Mer- 
rem, 1820) JX114025, Montivipera xanthine (Gray, 
1849) KJ415303 and Causus defilippi (Jan, 1863) 
AY223556). The final alignment comprised the 
910 bp cytb fragment. To choose the best model of 
molecular evolution (НКУ+1+С), we used Akaike's 
information criterion (AIC) in JMODELTEST 
2.1.1 (Darriba et al. 2012). Phylogenies were re- 
constructed using Bayesian inference (BI) and 
maximum likelihood (ML) approaches. BI for cytb 
data performed in MRBAYES 3.2.2 (Ronquist and 
Huelsenbeck 2003). Each BI analysis started with 
random trees and performed two independent runs 
with four Markov chains Monte Carlo (MCMC) for 
5 million generations with sampling every 1000th 
generation. Consensus trees constructed based on 
the trees sampled after the 2596 burn-in. ML analysis 
calculated with TREEFINDER (Jobb 2011). Boot- 
strap analysis employed 1000 replicates. Final trees 
obtained in FIGTREE v1.4.0 (http://tree.bio.ed.ac. 
uk/software/figtree/). 

The level of genetic differentiation in cytb based 
on p-distances was estimated in Mega 7 (Kumar et 
al. 2016). The final alignment for the calculation of 
p-distances comprised the 663 bp cytb fragment. 


RESULTS OF THE ANALYSES 


1. Morphological characters 


The taxonomical revisions of the shield-head vi- 
pers of the Great Caucasus (Orlov and Tuniyev 1986; 
Vedmederya et al. 1986; Nilson et al. 1994; Tuniyev 
and Ostrovskikh 2001) reveal the complicated pat- 
tern of cryptic speciation in two lineages of vipers — 
“kaznakovi”-complex and “ursinii”-complex, with 
morphologically hardly differentiated species mainly 
in the east sector of the Great Caucasus ( Tuniyev S. 
et al. 2012). A similar pattern is observed in the 
Transcaucasia, where discrete forms inhabit relative- 
ly limited territory, reflecting the prolonged isolation 
not only for historical reasons, but also due to human 
activity here since ancient times. On the plains and 
in the foothills of Eastern Georgia, the intensive use 
of land for agriculture and livestock has contributed 
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to the extreme fragmentation of habitats for animals 
whose survival is associated with the conservation of 
non-plowed steppes. 

The steppe vipers found in some localities of the 
Shirak Plateau in East Georgia seems to belong to 
relict micropopulations, which are analyzed below. 

A comprehensive analysis of external morphology, 
altitude gradient in habitats, typology of biotopes and 
genetic analysis revealed a high degree of similarity 
of populations of the steppe vipers from Azerbaijan, 
known as Pelias shemakhensis, and from East Geor- 
gia. These results allowed to consider newly discov- 
ered steppe vipers in Eastern Georgia within the taxa 
P. shemakhensis. However, according to the results of 
the comparison of the meristic characters between 
samples of P. shemakhensis from East Georgia and 
Azerbaijan the following differences were recorded 
(see Tables 3 and 4). 

Males. Significant differences were obtained for 
the number of subcaudals, which is higher in the 
vipers from Georgia. However, the overall trend is 
observed for pholidosis oligomerization from West to 
East from populations of East Georgia to the popula- 
tion from Azerbaijan. 

Females. The higher number of subcaudals in 
males was recorded in the vipers from Georgia too. 
In addition, the higher number of scales around mid- 
body (Sq.2) was noted in the vipers from Georgia. 


Table 3. Comparison of morphological characters of males of Pelias 
shemakhensis from populations from East Georgia and Azerbaijan. 
See Table 2 for abbreviations. 





Mean Mean 
P. shemakhensis P. shemakhensis 





Character East Georgia Azerbaijan t-value p 
3d 

Pr. 2.0 2.5 1.88982 20.05 
Ven. 135.2 134.0 1.22988 0.05 
9:6. 35.2 33.0 2.68373 <0.05 
Ар. 1.2 1.0 0.59761 >0.05 
C.s. 13.2 12:5 0.29111 — 20.05 
Sq2 214 21.0 0.59761 >0.05 
Supralab. 18.2 18.0 0.59761 >0.05 
Sublab. 19.8 20.0 0.59761 >0.05 
Ec. 19.8 17.5 1.82862 20.05 
ZZ 117.0 106.0 0.81582 20.05 
Lor. 9.8 9.0 0.52164 0.05 
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Table 4. Comparison of morphological characters of females of 
Pelias shemakhensis from populations from East Georgia and Azer- 
baijan. See Table 2 for abbreviations. 

















Mean Mean 
Character P. shemakhensis P. shemakhensis уййне 
East Georgia ^ Azerbaijan р 
99 PP 
Pr. 1.5 1.7 0.38069 >0.05 
Ven. 139.5 134.9 1.74815 20.05 
S.c. 29.0 26.3 2.34473 <0.05 
Ap. 1.0 0.9 0.43033 >0.05 
C.s. 16.0 13.3 112441 >0.05 
Sq2 21.5 21.0 2.88675 . «0.05 
Supralab. 19.5 18.0 8.66025 «0.001 
Sublab. 19.5 21.0 1.50756 — 20.05 
Ec. 18.5 18.8 0.23973 >0.05 
ZZ 109.5 109.6 0.01082 >0.05 
Lor. 13.5 10.1 210697 >0.05 
Tree Diagram for 4 Variables 
Single Linkage 
Euclidean distances 
P. eriwanensis 
P. lotievi east 
P. schemachensis 
Ge. 
Р. ^ 
Az. 
5 10 15 20 25 30 35 


Linkage Distance 


Fig. 2. Results of the cluster analysis (UPGMA method) of four 
groups of males of Pelias spp. based on pholidosis characters. 


The pholidosis oligomerization was registered from 
West to East: from populations of East Georgia to the 
population from Azerbaijan. 

With regard to the heterogeneity and the low 
number of individual samples, sexual dimorphism in 
the populations studied is expressed in different sets 
of characters. Below we consider these characters for 
each population (see Tables 5 and 6). 

Pelias shemakhensis, East Georgia (n-7). The 
number of subcaudal shields (S.c) and wings of zigzag 
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Tree Diagram for 4 Variables 
Single Linkage 
Euclidean distances 


P. eriwanensis 





P. lotievi east 





P. schemachensis 





10 12 14 16 18 20 2 24 
Linkage Distance 


Fig. 3. Results of the cluster analysis (UPGMA method) of four 
groups of females of the steppe vipers based on pholidosis char- 
acters. 


(ZZ) are higher in males than in females, while the 
number of ventral shields (Ventr.), supralabial shields 
(Supralab.) and number of dorsal scales around the 
posterior part of body (Sg3) is lower in males than 
in females. 

Pelias shemakhensis, Azerbaijan (n=12). The 
number of subcaudal shield is higher in males than 
in females. 

Due to differences in the method of counting of 
the number of crown shields (C.s.) and number of 
supralabials (Supralab.) by Kukushkin et al. (2012) 
from Azerbaijan and in our sample from Georgia, we 
eliminated methodologically different data about 
two indicated characters and utilized only eight mer- 
istic characters (Pr., Ven., S.c., Sg.2, Sublab., Ec., ZZ., 
Lor.) in the cluster and discriminant analyses. 

The cluster analysis based on eight meristic 
characters of pholidosis of the males has resulted in a 
dendrogram consisting of two big clusters: 1) vipers 
of mountain areas of the Transcaucasia and North 
Caucasus; and 2) vipers from the East-Transcaucasia 
Depression. The first cluster is further divided into 
subclusters of Pelias eriwanensis and P. lotievi, where- 
as the second one, into two subclusters, correspond- 
ing to two populations of P. shemakhensis (Fig. 2). 

The cluster analysis based on eight meristic char- 
acters of pholidosis of the females showed absolutely 
identical results with those of the males (Fig. 3). 

The geographical variability of morphological 
characters in populations of the steppe vipers was 
also analyzed by CDA, allowing a comparison of the 
preliminary selected groups in the complex of char- 
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Table 5. Sexual dimorphism of Pelias shemakhensis from East Georgia population. See Table 2 for abbreviations. 
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Can. 5 $51 2 2 0.59761 20.05 
Sal 5 ae 2 2 0.59761 >0.05 
Sq? 5 ie 2 sea 0.13894 >0.05 
"x Sa 
Supralab, 5 ги 2 йоз 3.04725 <0.05 
Sublab. 5 i 2 коо 0.70321 >0.05 
Be 2 oem 2 в 1.28593 е 
=з ою E әш ан 
Lor. 5 о 2 "our 2.52809 20.05 





acters (Tyurin et al. 2003). For comparison, Pelias 
lotievi east was taken as an outgroup. 

We use a complex of seven meristic characters 
(Pr, Ven., S.c., Sq.2, Sublab., Ec., Lor.) for which the 
reliable differences were obtained in the statistical 
analysis. The vipers were divided a priori into eight 
sexual and geographical groups. 

The results of CDA showed enough high accu- 
racy of the division of the geographical groups. The 
accuracy for the males is the following: Pelias eriwan- 
ensis — 77.8%, P. lotievi east — 100%, P. shemakhensis 
Ge — 100%, P. shemakhensis Az — 100%. 

All females showed lower accuracy of the division 
of geographical groups: Pelias eriwanensis — 46.2%, 


P. lotievi east — 73.726, P. shemakhensis Ge — 100%, 
P. shemakhensis Az — 90%. 

The results of CDA showed that in space of dis- 
criminant functions the males formed three groups 
(Fig. 4): 1) Pelias eriwanensis from Armenia; 2) P. sh- 
emakhensis from Azerbaijan and Georgia; 3) P. lotievi 
east from Dagestan. The separation ofthe species level 
on the first discriminant function is present between 
P. lotievi east and P. shemakhensis of Azerbaijan and 
Georgia, while division P. lotievi east and P. eriwanen- 
sis is more clear for the second discriminant function. 

The distribution of the females in space of 
discriminant functions (Fig. 5) looks more hetero- 
geneous with the formation of three independent 
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Fig. 4.  Two-dimensional 
scatterplot of samples of 
males of the steppe vipers in 
space of CDA function on 
the complex of morphometric 
characters. 


Fig. 5. Two-dimensional 
scatterplot of samples of 
females of the steppe vipers 
in space of CDA function on 
the complex of morphometric 
characters. 
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Table 6. Sexual dimorphism of Pelias shemakhensis from Azerbaijan population. See Table 2 for abbreviations. 
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Xtm Xtm 
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Table 7. Mahalonobis distances (right) and significance levels (left) between samples of males of Pelias eriwanensis, P. lotievi east, P. shem- 
akhensis Ge. and P. shemakhensis Az., according to CDA results. 








P. eriwanensis P. lotievi east P. shemakhensis Ge. P. shemakhensis Az. 
P. eriwanensis 10.36361 12.90747 22.24767 
P. lotievi east 0.000718 26.58053 41.39274 
P. shemakhensis Ge. 0.005897 0.000018 6.06630 
P. shemakhensis Az. 0.011396 0.000181 0.599716 





Table 8. Mahalonobis distances (right) and significance levels (left) between samples of females of Pelias eriwanensis, P. lotievi east, P. she- 
makhensis Ge. and P. shemakhensis Az., according to CDA results. 





P. eriwanensis P. lotievi east P. shemakhensis Ge. P. shemakhensis Az. 
P. eriwanensis 3.05010 11.78868 7.679577 
P. lotievi east 0.034834 15.79790 8.154632 
P. shemakhensis Ge. 0.063467 0.013326 9.672418 


P. shemakhensis Az. 0.000944 0.000187 0.145138 
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Table 9. Contribution of different morphological characters in 
separation of males of Pelias eriwanensis, P. lotievi east, P. shem- 
akhensis Ge., P. shemakhensis Az., according to CDA results. See 
Table 2 for abbreviations. 








Character Standardized coefficient of Character 
the first discriminant function rank 

Pr. 0.25 4 
Ven. 0.14 1 
Sc. 0.33 3 
Sq.2 0.25 5 
Sublab. 0.45 2 
Fc. 0.23 6 
ZZ 1.03 1 
Lor. 0.10 8 





groups. At that Pelias shemakhensis from Azerbaijan 
and P. shemakhensis from Georgia are well separated 
in both discriminant functions, whereas P. lotievi east 
and P. eriwanensis are separated only in the second 
discriminant function. 

The species level differences (on the first dis- 
criminant function) are noted in Pelias lotievi and 
P. shemakhensis (combined group from Georgia 
and Azerbaijan), whereas P. lotievi east and P. eri- 
wanensis are characterized rather by the subspecies 
level differences (more clear in the females). Finally, 
P. shemakhensis from Azerbaijan and P. shemakhensis 
from Georgia demonstrate only weak population 
level differences, perhaps due to the small number of 
samples. 

The obtained results confirm the high degree of 
morphological separation of the compared samples of 
vipers. Degree of likeness between the selected sam- 
ples in CDA was estimated as the value of the Makha- 
lonobis distance (Tyurin et al. 2003). The distances 
between the centers of the samples of the males of vi- 
pers varied from 6.1 to 41.4. The minimum value was 
shown between the males of Pelias shemakhensis from 
Georgia and P. shemakhensis from Azerbaijan (6.1) 
and between P. lotievi east and P. eriwanensis (10.4), 
whereas the maximal value (41.4) — between the 
males of P. shemakhensis from Azerbaijan and P. lotievi 
east from Dagestan (Table 7). For the females, this dis- 
tance between the centers of the samples varied from 
3.1 to 15.8. The minimum (3.1) distance was recorded 
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Table 10. Contribution of different morphological charac- 
ters in separation of females of Pelias eriwanensis, P. lotievi 
east, P. shemakhensis Ge., P. shemakhensis Az., according to 
CDA results. See Table 2 for abbreviations. 








Character Standardized coefficient of Character 
the first discriminant function rank 

Pr. 0.49 2 
Ven. 0.01 8 
S.c. 0.06 7 
Sq.2 0.08 6 
Sublab. 0.45 3 
Fe. 0.22 5 
22 0.91 1 
Lor. 0.24 4 





between the females of Р. lotievi east и Р. eriwanensis 
and P. shemakhensis from Azerbaijan and P. eriwan- 
ensis (8.2); the maximal distance (15.8) — between 
females P. shemakhensis from Azerbaijan and P. lotievi 
east from Dagestan (Table 8). 

The contribution of different morphological char- 
acters to discrimination of the groups is different. 
Because the first discriminant function takes into 
account the most percent of the dispersion, and the 
division of the vipers into main groups occurs exactly 
along it, the contribution of characters to the divi- 
sion of the groups is described based on values of this 
function (Tables 9 and 10). 

The followings characters made a maximal contri- 
bution to discrimination of the groups of the males 
(Table 9): wings of zigzag (ZZ), sublabials, and sub- 
caudals. The followings characters made a maximal 
contribution to discrimination of the groups of the 
females (Table 10): wings of zigzag (ZZ), preventral 
shields (Pr.), and sublabials. 

As a result of our study, new information was 
obtained about the morphological characters and 
geographical variability of the three viper species. 
Substantial differences in the mean values in a num- 
ber of metric and meristic characters of the vipers are 
most valuable (Table 3 and 4). The discrimination of 
the three groups into four samples, selected on the 
principle of geographical and sexual identity using a 
cluster analysis and CDA (Figs. 2—5; Tables 6-10) is 
also significant. 
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Table 11. Uncorrected p-distances between samples of Pelias spp. 


99 








P. lotievi P. renardi P. ebneri P. shemakhensis P. eriwanensis 
P. lotievi 0.006 
P. renardi 0.011 0.005 
P. ebneri 0.023 0.021 0.006 
P. shemakhensis 0.025 0.022 0.017 0.002 
P. eriwanensis 0.028 0.025 0.020 0.011 0.004 





2. Genetic analysis 


Cytochrome b tree (Fig. 6) shows the following 
structure. The main cluster includes Pelias lotievi, 
P.renardi Christoph, 1861 and highly supported clus- 
ter (0.96 BI and 80.3 ML) that unites three species: 
P. ebneri (0.84 BI and 71.4 ML), P. shemakhensis (1.0 
BI and 99.7 ML) and P. eriwanensis (1.0 BI and 97.5 
ML). The P. eriwanensis cluster is subdivided into 
subclusters according to the local populations. 

The p-distance between Pelias shemakhensis and 
P. lotievi is 2.5%, and between P. shemakhensis and 
P. eriwanensis is 1.1 % (Table 11). 

The results allow to conclude that the taxonomic 
differences vary in all four groups, three of which 
(Pelias eriwanensis, P. lotievi east, P. shemakhen- 
sis) differ on the species level, and two of them (P. 
shemakhensis Ge, P. shemakhensis Az), on the subspe- 
cies level. 


SYSTEMATICS 
Family Viperidae Laurenti, 1768 
Genus Pelias Merrem, 1820 


Pelias shemakhensis kakhetiensis Tuniyev, 
Iremashvili, Petrova et Kravchenko ssp. nov. 


Holotype. SNP 1059, vicinity of the Dedoplisz- 
karo Village, Fortress Khornabudji, Shirak Plateau 
(Kakhetia, East Georgia), 06.2016, Coll.: N. Seturid- 
ze (Fig. 7). 

Paratypes. Six specimens: SNP 1052 (4 juv.), 
vicinity of the Chinkani Village, Shirak Plateau 
(Kakhetia, East Georgia), 09.2015, Coll.: G.N. Ire- 
mashvili (Figs. 5 and 6); SNP 1053 (2 specimens, 
adult mail and adult female), vicinity of the Chinkani 
Village, Shirak Plateau (Kakhetia, East Georgia), 
04.2016, Coll.: G.N. Iremashvili (Fig. 8A, B). 


Diagnosis. Small-sized snake, different from 
the nominative subspecies by larger numbers of 
pholidosis characters, first of all, maximum number 
of subcaudals and supralabials. Apical shields are 
developed in all specimens, in contrast to the nomi- 
native subspecies, some specimens of which have 
a combination of four canthals in absence of apical 
(a character unique for the whole Pelias complex). 
A weakly undulated brown zigzag is located in a gray 
or brown dorsal background. In adults, extra-large 
patches and some minor dots on the sides of the body 
can be present in the posterior third of the body, 
unlike the nominative subspecies with small dots in 
2-3 staggered rows along the whole body (Fig. 9). 
Crown spot is triangular in shape with a contoured 
white stripe in front of the frontal shield. Adults have 
a white stripe becoming less contrast, with the con- 
tinuing triangular crown spot. 

Description of the holotype. Adult male having 
the following morphological characters: total length 
(Т1) – 350 mm, length of tail (L.cd.) — 35 mm. 
Number of preventrals — 2; ventrals — 134; subcau- 
dals — 35; apical — 1; canthals — 5; crown shields 
(C.s.) — 5. Upper preocular shield is in contact with 
nasal. Number of rows of dorsal scales around a neck 
and at midbody — by 21, and at posterior part of the 
body — 17, sublabials on the right — 10, on the left — 9, 
supralabial shields on both sides — 9, shields around 
eyes — 10 from the right and 8 from the left, wings of 
zigzag — 114, loreals — 5 from the right and 4 from the 
left. Length of the head — 20.3 mm, width — 12.6 mm, 
height — 8.9 mm, length of pileus (Pil.) — 12.6 mm. 
Frontal, nasal and parietal shields are not divided. 
The holotype specimen is grey from above with dark- 
brown zigzag, whose wings are rounded. On the sides 
of the body there are extra-large patches, whereas 
several minor dots are present in the posterior third 
of the body. Pattern of the head not merges with the 
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144 P. lotievi (9) 
160 Р lotievi (8) 
422 P. lotievi (12) 
243 P. lotievi (10) 
208 P. renardi (11) 
186 P. renardi (14) 
187 P. renardi (14) 
185 P. renardi (14) 
192 P. renardi (14) 
201 P. renardi (13) 
204 P. renardi (13) 
P. ebneri KC176745 
P. ebneri FR727093 
P. ebneri FR727094 
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402 P. shemakhensis (7) 
404 P. shemakhensis (7) 
399 P. shemakhensis (7) 
420 P. shemakhensis (6) 
171 P. eriwanensis (1) 


















172 P. eriwanensis (1) 
170 P. eriwanensis (1) 
42 P. eriwanensis (4) 
219 P. eriwanensis (4) 
220 P. eriwanensis (3) 
43 P. eriwanensis (2) 

45 P. eriwanensis (5) 





L Nc OUTGROUP 
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Fig. 6. BI tree for the Pelias spp constructed using cytb haplotypes. Bootstrap values are as follows: Bayesian posterior probabilities / 
maximun likelihood bootstrap supports. For specimen codes, see Table 1. For location numbers (in brackets), see Fig. 1. 


New shield-head viper from Georgia 





ДЫ бе". о РА 
Fig. 7. Holotype of Pelias shemakhensis kakhetiensis ssp. nov. (SNP 
1059). 
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Fig. 8. Part of the type series of Pelias shemakhensis kakhetiensis 
ssp. nov. (SNP 1053), left — male, right — female: A — view from 
above, B — view from below. 


zigzag. The crown spot is triangular, contoured by 
white stripe in front of the frontal shield. Coloring of 
the throat and lower jaw is light, almost white. Belly 
and low part of the tail are dark, almost black: the tip 
of the tail is yellow. 

Description of the paratypes. The paratypes are 
corresponding to description of the holotype with in- 
significant variations in a size and meristic characters 
(Table 12). 

The main background of adult males is grey, fe- 
males — grey or pinkish-grey. Zigzag in both sexes 
continuous, undulating dark brown. Zigzag is sepa- 
rated from pattern of the head. Coloration of juve- 
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Fig. 9. Nominative subspecies of Pelias shemakhensis (Azerbaijan, 
vicinity of the Maraza Village) with typical rows of small dots on 
lateral sides of body. 





Fig. 10. Color and patterns of the juvenile specimen of Pelias shem- 
akhensis kakhetiensis ssp. nov. (paratype, SNP 1052). 


niles is identical to that of the adults, but the char- 
acteristic triangular dark patch in the anterior part 
of the frontal shield with front white stripe, outlined 
more contrast (Figs. 10 and 11). The coloration of the 
pileus is dark-brown, almost black in the males. On 
the sides of the body of the juveniles, there are dots 
in three rows or one row of dark spots. The adults 
have extra-large patches on the sides of the body with 
the second row of several minor dots presented in the 
posterior third of the body. Coloration of the throat 
and lower jaw is light, almost white in the adults. The 
belly in young and adult females is dark spotted, with 
a predominance of dark tones, whereas in males of 
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Fig. 11. Heads of two paratype specimens of Pelias shemakhensis 
kakhetiensis ssp. nov. (SNP 1052; A and B) showing a characteris- 
tic triangular dark patch in the anterior part of the frontal shield 
defined in front by white stripe. 





Fig. 12. Distribution of Pelias shemakhensis kakhetiensis ssp. nov. 
(1—5) and P. shemakhensis shemakhensis (6—11): 1 — Norio Village; 
2 — Artsiviskheoba Mountain; 3 – Khornabudji Fortress (Quin 
Tamar); 4 — Chinkani Village; 5 — Vicinity of the Zemo Kedi 
Village; 6 — Karatala Village near Nukhi (= Sheki) Town; 7 — 
Demirchi Village; 8 — Khynysly Village; 9 — Angekharan Village; 
10 — Shemakha Town, 11 — Maraza Village. 


all ages it is almost black; the tip of the tail in both 
sexes and ages is yellow from below. Dark strips are 
present, or not in pattern of the front edge of the su- 
pralabials; sublabial shields have pattern, reminiscent 
of the “temminck” figure, but not as contrast as those 
of Pelias renardi. 

Etymology. The subspecies name is after histori- 
cal region of Georgia — Kakheti, where a large part of 
the range is located. 

Geographical distribution and biotopes. The 
distribution range of the species in Georgia (Fig. 12) 
covers the peripheral areas of the Shirak Plateau in 
two neighboring regions — vicinity of the Chinkani 
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Fig. 13. Killed specimen of Pelias shemakhensis kakhetiensis ssp. 
nov., Artsivis Kheoba Mountain. 


Village and Zemo Kedi Village; in vicinity of the De- 
dopliszkaro Village near the Khornabudji Fortress, 
also known as the Fortress of Queen Tamar; at north 
foot of the mountain Artsivis Kheoba (Fig. 13). In 
the collection of the SMG, there is specimen from the 
vicinity of the Norio Village (SMG 125, year 1933, 
unnamed collector), which represents the western 
limit of the distribution of the species in Georgia 
(Fig. 14). 

Recent distribution of the nominative subspecies 
is limited by the Shemakha District of Azerbaijan, 
from where it is known from vicinity of Shemakha 
Town, Demirchi Village; Gobustan District — Mara- 
za Village. Alekperov (1978) noted records of the 
nominative subspecies in the town of Shemakha, be- 
tween Shemakha and Angekharan Village, between 
Shemakha and Khynysly Village, at the foot of the 
Pirdirechidag Mountain, on the slopes of the Gyzka- 
lasy Mountain. For the first time, it was recorded in 
Azerbaijan near Karatala Village in vicinity of town 
Nukhi (Nikolsky, 1913), formerly a part of Georgia. 

Both preserved fragments of its former distribu- 
tion area in Georgia and Azerbaijan are clearly the 
relic (Fig. 12). 

Like in Azerbaijan, the biotopes of the species are 
represented by shibliaks of Paliurus spina-christi, eco- 
tones of dry Andropogon steppes and derivates of the 
broad-leaved forests on the stony slopes of foothills 
of the Shirak Plateau in East Georgia (Fig. 15), since 
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the entire central part of the Shirak Plateau is opened 
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Fig. 14. Pelias shemakhensis kakhetiensis ssp. nov. (SMG 125) from 
vicinity of the Norio Village. 





Fig. 15. Biotopes of Pelias shemakhensis kakhetiensis ssp. nov. near 
the Chinkany Village: shibliaks of Paliurus spinachristy and eco- 
tones of dry Andropogon steppes and derivates of the broad-leaved 
forests. 
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Fig. 16. The entire central part of the Shirak Plateau (Didi Shi- 
raki) plowed and cultivated for crops. 


solute minimum is —22°C, with an average absolute 
minimum of —14°C. The average annual rainfall is 
450—600 mm. 

The Shirak Plateau has an altitude of 865 m above 
51. with the predominance of wave-hill terrain and 
steep dumping of the Alazani Valley. The average 
temperature in January is 3.2°C, with the absolute 
minimum of —32”C. The average temperature in July 
is 22.8°C, the absolute maximum is 38°C. The average 
annual temperature is 10.3°C. The average annual 
precipitation is 490 mm. 

Both plateaus are currently woodless, derivatives 
of forests are found on the steep slopes of the bound- 
ary areas, whereas the plain parts of the plateaus are 
covered by shibliaks of Paliurus spina-christy, second- 
ary dry Andropogon steppes and agricultural crops. 

The Shirak Plateau is separated from the Shem- 
akha-Gjurdzhivan Plateau by the Alazani-Agrichaj 
Plain, with a length of about 300 km, with a width 
of 20—25 km. In the past, it was fully covered with 
deciduous forests, which was a natural barrier for 
penetration of the steppe and semi-desert species of 
plants and animals. 

In the Shirak Plateau, Pelias shemakhensis prefers 
to inhabit the most elevated northern peripheral 
sites, the fact also noted by E.S. Balavadze (1966). 
This author indicated that the steppe viper belongs 
to the reptiles, inhabiting the specific zone and not 
penetrating to another sites on the Shirak Plateau. 
Balavadze (1966) called the steppe viper as “ecologi- 
cal endemic”. According to our observations, the spe- 
cies is also found on the rocky slopes of the northern 
exposure, to the Alazani Valley. 
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As was reported previously (Tuniyev S. et al. 
2013), on the Shemakha-Gjurdzhivan Plateau, the 
species also inhabits the rocky slopes covered by 
shibliaks of Paliurus spina-christi and secondary dry 
Andropogon steppe. 

In biology of Pelias shemakhensis, we marked a 
unigue feature missing from other representatives of 
the “ursinii”-complex. In summer, the animals change 
to crepuscular-night activity, which is associated 
with extremely high daytime air temperatures ex- 
ceeding 35°C. So, that adult female (not included in 
the type series) was observed at 0—30 a.m. (Fig. 17). 
Alekperov (1978) noted twilight activity of these 
vipers in the vicinity of Shemakha. 

The statistical analysis of the external morpho- 
logical pholidosis characters shows relatively small 
differences between animals of two geographical 
clusters, which are largely determined by the small 
samples. However, long enough geographic isolation 
has resulted in the formation of the above marked 
differences, including elements of the coloring of the 
body and the head, which allows considering the vi- 
pers from the Shirak Plateau in Georgia and the She- 
makha-Gjurdzhivan Plateau in Azerbaijan as distinct 
subspecies. So, Pelias shemakhensis kakhetiensis ssp. 
nov. does not demonstrate a unique combination of 
apical and canthal shields missing in other members 
of the “ursinii complex”, but typical for specimens of 
the nominative subspecies. P. shemakhensis kakhe- 
tiensis ssp. nov. has a characteristic triangular dark 
patch in the anterior part of the frontal shield with 
a white stripe at the front, whereas the nominative 
subspecies has a fuzzy dark spot located in the center 
of the frontal shield. The white stripe in the anterior 
part of this shield is missing. Differences in the col- 
oration of the body and pholidosis are noted above. 

Probably, in the Pliocene, Pelias shemakhensis had 
a single area, stretching in the steppe and forest-steppe 
landscapes north of the Transcaucasian Depression 
from the vicinity of Tbilisi to southeast of the Eastern 
Caucasus foots. It is known that in the Paleogene of 
the Caucasus already occur all recent families of the 
reptiles (Bakradze and Chkhikvadze 1988); thus on 
the Caucasus and in the nearest environs of this Isth- 
mus lived both mesophylous and various xerophylous 
species of herpetofauna, identical or closely related 
to recent species of the Caucasus. Darevsky (1963) 
considered that by the end of the Pliocene in the Cau- 
casus the primary core of its recent herpetofauna had 
been formed with such genera as Agama (= aralau- 
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Fig. 17. Adult female of Pelias shemakhensis kakhetiensis ssp. nov. 
observed at 0—30 a.m. near Chinkani Village. 


dakia), Lacerta sensu lato (Darevskia, Lacerta sensu 
stricto), Ophisaurus (= Pseudopus), Anguis, Typhlops, 
Malpolon, and Vipera. 

The early orogenesis processes intensified in the 
late Neogene-Quaternary resulted in upraising of 
flattened surfaces and their dismemberment by can- 
yons. In this period the Pelias shemakhensis range was 
divided by the Alazani-Agrichaj Plain, completely 
covered by forest. The further fragmentation of 
habitats in both clusters took place solely under the 
influence of human activities, leading to exceptional 
recent rarities, both the nominative subspecies, and 
the newly described Pelias shemakhensis kakhetiensis 
ssp. nov. The inclusion in the Red List of IUCN and 
Red Data Book of Azerbaijan with status category 
CR Blab (i, ii, iii, iv, v) has been requested for the 
nominative subspecies (Tuniyev S. et al. 2013). The 
similar status for Pelias shemakhensis kakhetiensis ssp. 
nov. is required in Georgia. 

The main intrigue around Pelias shemakhensis 
remains in its origin and relationship with other 
members of the “ursinii”-complex in the Caucasian 
Isthmus. These questions can be answered based on 
the results of the cluster, discriminant and genetic 
analyses provided here. 

For the description of Pelias shemakhensis, we 
used only museum specimens from Azerbaijan (Ku- 
kushkin et al. 2012) without genetic analysis. 

Based on the results of the cluster analysis of 
pholidosis characters (Tuniyev S. et al. 2013), the 
minimum distance was noted between clades of Pelias 
renardi east, P. renardi west and P. lotievi, whereas 
P. shemakhensis and P. eriwanensis form sister clus- 
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ters. However, results of the discriminant analysis of 
the evolutionary conservative females suggest that 
P. shemakhensis (Azerbaijan) and P. lotievi have sub- 
species rather than species level differences, whereas 
P. eriwanensis and P. shemakhensis are different on the 
species level (Tuniyev S. et al. 2013). 

Inclusion of the materials from Georgia in the 
cluster, discriminant and genetic analyses in this 
paper changed the picture. Based on the results of 
the cluster and discriminant analyses, Pelias eriwa- 
nensis and P. lotievi should be given a subspecies rank, 
whereas P. shemakhensis clearly deserves a species 
rank. Results of the genetic analysis are opposite: 
P. shemakhensis and P. eriwanensis are considered as 
sister species. 

At present, it is difficult to talk about relation- 
ships and distribution routes of ancestors of the 
entire complex of the steppe vipers. Several conflict- 
ing hypotheses were suggested (Zinenko et al. 2011; 
Ferchaud et al. 2012; Gvozdik et al. 2012; Tuniyev S. 
et al. 2013) about the routes of distribution of the 
steppe vipers. 

We continue to share the idea about autochtho- 
nous origin for small shield-head vipers of the Cau- 
casus ( Tuniyev et al. 2013) taking into account their 
ancient origin in the Caucasian Ecoregion and as- 
tonished variety of forms of the “kaznakovi”-complex 
and "ursinii"-complex both on the Great Caucasus 
and in the Transcaucasia. 

Nilson and Andrén (2001) are very close to this 
idea, suggesting the hypothetical scenario of origin 
and distribution of the steppe vipers from the north- 
east to the southwest and considering the Caucasus 
as one of Pliocene refugia for the members of the 
“ursinii” sensu lato. Possibly, genetic analysis of the 
nuclear DNA will allow coming closer to understand- 
ing this issue. 
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